RATE OF DECOHERENCE FOR AN ELECTRON 
WEAKLY COUPLED TO A PHONON GAS 



RICCARDO ADAMI AND LASZLO ERDOS 

Abstract. We study the dynamics of an electron weakly coupled to a phonon gas. The initial state of the 
electron is the superposition of two spatially localized distant bumps moving towards each other, and the 
phonons are in a thermal state. We investigate the dynamics of the system in the kinetic regime and show 
that the time evolution makes the non-diagonal terms of the density matrix of the electron decay, destroying 
the interference between the two bumps. We show that such a damping effect is exponential in time, and the 
related decay rate is proportional to the total scattering cross section of the electron-phonon interaction. 



1. Introduction and result 

Quantum interference of a state can be destroyed via interactions with an external environment. 
This phenomenon is called decoherence and it is recognized to be responsible for the transition from 
quantum superpositions to classical probability measures. This feature makes it relevant in the field of 
quantum information, where the superposition of states allows for parallel computation and therefore 
it must be preserved. 

Rigorous results on decoherence have so far concerned systems with a discrete energy spectrum 
or decoherence induced by practically instantaneous scattering at time zero. For the first model, we 
mention that the decoherence of an electron trapped in a harmonic oscillator and interacting with an 
electromagnetic field in the dipole approximation was studied in |DS| . For the second situation, in 
|DFT| IAFFTl| [AFFT2J it was showed that decoherence is the most relevant effect in the dynamics 
of many heavy particles in an environment of light particles, when the mass ratio is small. For an 
extensive physical introduction to decoherence we refer to jGJKKSZ, ID]. 

Here we investigate the occurrence of decoherence in the dynamics of an electron weakly coupled 
to a phonon field. The energy spectrum of the electron is continuous and the phonons constantly 
bombard the electron along its time evolution. We will study the problem in the kinetic regime with 
a weak coupling. We introduce a scale parameter e and analyze the time evolution at large space and 
time scales of order e _1 . If (x, t) denote the original (microscopic) space and time coordinates, then 
we introduce the macroscopic space and time coordinates (X, T) by 

T := et, X := ex. 

The coupling constant between the electron and phonon field will be given by A = \fe. 

On the macroscopic space and time scales, the dynamics of the electron in such a limit is described by 
a linear Boltzmann equation, as proven in [Ej. In this paper we prove that, in addition to such behavior, 
on the microscopic scale the quantum interference remains detectable and decays exponentially in time 
with a rate growing with the strength of the electron-phonon interaction. 
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To be more specific, we choose a special initial state that exhibits a strong interference pattern 
without phonon interaction. We assume that initially the electron lies in the superposition of two 
distant macroscopic wave packets ("bumps") in R d (d > 3) moving towards each other, namely 

^g(x) := ipo,+ (x) + Vo,-( x ) 

(1.1) ^ >+ (*):= e d l 2 f{ex + Qy Px 

^_{x) := e d / 2 f(ex-Q)e- lP - x 

where / is L 2 -normalized and satisfies the following hypothesis of regularity 

(1.2) / G H% +6 {R d ). 

The vectors P, Q e R d are parallel, i.e. P-Q = \P\\Q\ and assume The wave packets 4>q _ and 

tpQ + are localized around the macroscopic points Q and — Q and they move towards each other by 
a momentum P. The function / describes the macroscopic envelope of the bumps. We remark that 
ipQ is not normalized in L 2 , but it is easily seen that lim e ^o HV'oll = v^j so we wm & l wa ys consider e 
sufficiently small to guarantee 1 < HV'olk < 2. 

To explain our result, we preliminarily examine what happens if the interaction between electron and 
field is absent. In this case the state of the electron ipf = ipf ivee evolves according to the free Schrodinger 
equation idtipt = — |^^t- The two wave packets maximally overlap at time t = e^lQI/l-Pl, when 
the probability density of finding the electron at x, apart from the normalization factor, is given by 



\iPI Me (x)\ 2 = 2{1 + C os{2P ■ x))\e d ' 2 (e^ f) (ex) 



The shape of the bumps has undergone a semiclassical change and their supports are shifted to the 
origin. If the state were classical, the probability densities of the bumps should add up, yielding the 
factor of 2. The cosine term carries the typical interference fringes and therefore the information that 
the electron lies in a quantum superposition of two wave packets. 

Let us remark the presence of two spatial scales in l^ff e (^)| 2 : the envelope of the wave function 
is spread on a length of order while the interference pattern manifests itself on a length scale 
of order 1. The size of the interference fringes can be quantified evaluating the Fourier transform of 
1^1 free 1 2 an< ^ then performing the limit e — > 0. One easily sees 

. _ ( 2 ifP = 

(1.3) lim / d*e 2 ^|Vf >free (*)| 2 =1 1 if P = ±P 

^ otherwise. 

If, instead of being isolated, the electron interacts from the beginning with a phonon field that 
initially lies in a thermal state, then the interference fringes are expected to be damped by a factor 
exponentially decreasing in the collision time in macroscopic units, i.e. with |Q|/|P|. 

Indeed, we prove that 

CTplQI 

^±2iP-XF( ■ — 



lim J dxe ±m pf{x) = e~^T 

where ap, given in (|7.4p . is the total scattering cross section for an electron with velocity P interacting 
with the phonon field, and pf(x) is the probability density of finding the interacting electron at x at 
time t. The interacting Hamiltonian will be defined in Section HI Our method easily extends to more 
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general initial states whose free evolution exhibits interference fringes. For simplicity we discuss only 
the particular case given by (jl.ll) . 

To state our result in a more precise way, we use the Wigner function formalism instead of the wave- 
function description of the state of the electron. At time zero, the Wigner function of the electron 
reads 

dy 



(1.4) W (x,v) := / —^- e - w ^l{x + y/2)%{x-y/2). 

According to the decomposition of ipQ given in (|1.1|) . formula (|1.4|) yields a decomposition in four terms 
for the Wigner function of the initial state 

W = HW, W ,aa'(x,v) := f e- iv -y^ a (x + y/2)r 0>a ,(x - y/2) , 

where the so-called diagonal terms Wo,++ and TVo, represent an electron in the state V'o + an d V'o - 

respectively, and the non-diagonal terms Wq^ and Wo t |_ represent the interference between the 

states ipQ , and ipQ _ . 

Since the time evolution of the Wigner function is linear, one can study separately the evolution of 

diagonal terms W ++ {t) and W (t) and of non-diagonal ones (i) and W (_(t). We will use the 

notations W± i ±j ree (t) and W±±(t) for the corresponding components of the Wigner transform of the 
state at time t under the free and the interacting evolutions, respectively. The interacting evolution 
will be defined in Section [H We remark that after tracing out the phonon degrees of freedom, the 
electron will be in a mixed state and thus W±±(t) will be Wigner functions of density matrices. 

In order to easily define the phonon field operators in a rigorous way we consider the system confined 
in a d-dimensional box A^ = [— L/2, L/2] d and then perform the thermodynamic limit L — > oo. Such a 
limit is taken after the phonon trace but before the scaling limit e — > 0, and the result is uniform in L. 
The parameter L will be suppressed in the notation of W±,±(i). Before performing the thermodynamic 
limit L — » oo, any integral in the space variables will be thought of as one over the volume A, and any 
integral in the momentum variable is actually a summation over A* = (L _1 Z) d . For the computation 
of the thermodynamic limit we refer to section 4.5 in [Ej. The free evolution will always be considered 
in the infinite volume to avoid taking the unnecessary thermodynamic limit in W±±f iee (t). 

We will test the evolved Wigner functions against observables, J(x,v) = J £ (x,v) that may scale 
with e to detect the interference fringes. Let (J,W) = f J(x,v)W(x,v)dxdv denote the expectation 
value of the observable J on a state given by the Wigner function W. 

We will always assume that 



(1.5) ||J||:=sup / sup | J e (£, v)\d£ < oo , 

e J v 

where we use the convention that the hat on functions defined on the electronic phase space denotes 
the Fourier transform in the space variable only, i.e. 



In the following Theorem ll.il we compare the kinetic limit for W+_(i) and W_ + (t) with the corre- 
sponding terms under the free evolution. With these notations, we have the following main theorem: 
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Theorem 1.1. Let the initial state of the electron be given by {Z3P with and fix a macroscopic 

time T . At any time t, let W-\ free(i) be a non-diagonal component of the Wigner function of a 

non-interacting ("isolated") electron evolved under the free evolution. Let W-\ it) be the analogous 

non-diagonal component of the Wigner function describing the electron at time t interacting with a 
phonon field through the interaction Hamiltonian ( f^.i[ ). Then, for any observable J £ satisfying (11. 5ft . 
we have 



(1.6) lim 



lim {J e ,W+-(£- 1 T)) - e- T ^(J £ ,W + _ ifree ( £ - 1 T)) 

L— >oo 



where ap, given in (|7.4p . is the total cross section for an electron with momentum P in the phonon 
field. 

To detect the destruction of the interference fringes, the two-scale structure of (t) and W ^(t) 

obliges us to test (t) and W |_ (t) against observables endowed with the same two-scale structure. 

A possible class of such observables J £ will be given in Section T4.41 In particular, for these observables 

one easily obtains that lim e _*.o(^ej W-\ ,frcc(£ -1 T)) exists. The precise statement is formulated in 

Corollary 14.11 

The evolution of the diagonal terms, W++ and W , on the macroscopic scales is given by the 

linear Boltzmann equation, 

(1.7) (d T + Ve(V) ■V x )Ft(X,V) = J a{V,U)F T (X,U)dU - ( J a{U,V)dU^jF T (X,V) 

with collision kernel o~(U,V) and free dispersion relation e(V) = \V 2 . The two terms on the right 
hand side are called the gain and loss terms, respectively. For the precise statement, see |Ej. 

Theorem 11.11 describes the evolution of the off-diagonal terms. The physical scenario can be ex- 
plained as follows: in the kinetic limit, the two initial bumps evolve independently as classical phase 
space probability densities obeying the linear Boltzmann equation (jl.7p . In addition to that, there is 
an interference term that equals the corresponding term for the free case, except for a damping factor, 
exponential in time. In other words, the interference term evolves according to a Boltzmann equation 
without gain term. 

From the technical point of view, we expand the time evolution into a Duhamel sum, trace out the 
degrees of freedom of the phonon and lastly perform the scaling limit. This is conveniently made using 
Feynman graph expansions, similarly to [EY] and [Ej. The novelty in our model is that the relevant 
graphs are no longer given by the ladder terms: these ones vanish in the limit and the resulting 
dynamics is given by the renormalized free evolution. 

The paper is organized as follows. In Section [2] we fix some basic notation and conventions; in 
Section [3] we state regularity assumptions on the dispersion relations and on the form factor of the 
interaction, and give a lemma on some estimates that will be used repeatedly along the paper. In 
Section [4] we define the model and the scaling, and specify the assumptions on the initial state and 
the observables. In Section [5] we explain the link with the result provided in [EJ. In Section [6] we 
prove that the contribution of the ladder graphs is negligible in the kinetic limit, and in Section [7] 
we compute the contribution due to renormalized free propagators. The paper ends by an appendix 
containing the proof of the lemma stated in Section (3) 
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2. Notation and conventions 

We model the electron as a spinless particle, so the state space for the electron is TL e '■= Lp er (A.£,). 
The electron dispersion relation is denoted by e(k), k £ and the Hamiltonian H e of a free electron 
acts on 7i e as follows 

(2.1) H c = e(-*V). 

One can think of the classical dispersion relation e(p) = j^j or the pseudo-relativistic one e(p) = 
\/\p\ 2 + M 2 , where M is the mass of the electron. 

Let 7 be the density operator representing the state of the electron; 7 is a positive operator on 
L 2 (M. d ) and its operator kernel is denoted by j(x,y). Let 7 be the Fourier transform of 7 as an 
operator, i.e. its kernel is given by 

7(p,«) := / ^e-^ + ^ 7 (x,y). 

The Wigner transform of 7 is defined as 

(2.2) W l(x ,v) ■*.,(,+ !,, _|) J*,. 

W 7 is a real function defined on the classical phase space of the electron. We recall the convention 
that the hat on functions defined on the phase space of the electron denotes Fourier transform in the 
space variable only, i.e. 

(2.3) W^v) = J J^ e -^W y (x,v) = (2*)"* 7 (v + \,v - |) . 

The rescaled (macroscopic) Wigner function is defined as follows 

' X 

e 

If the density matrix 7 is the orthogonal projection on the space spanned by the funtion ip (i.e. the 
electron lies in a pure state), then we denote the associated Wigner function by W v . We will sometimes 
use symbols like W(t), without explicitly referring to a density matrix 7. In such cases, we will always 
mean the Wigner function of the electron at time t, i.e. the Wigner transform of the reduced density 
matrix for the electron 7 ei to be defined in (|4.6p . 



W £ JX,V) = e~ d W^ ( 9,v 



The pure states of the phonon field are represented by vectors in the bosonic Fock space TC p h = 
[i 2 (A)] (8 " n , where ® s n is the n-fold symmetrized tensor product. We introduce the phonon 
creation and annihilation operators, at, a^, with momentum k, satisfying the usual commutation 
relations 

[a k ,a\,} = S(k - k') . 
The number operator of the phonons in mode k is 

J\f k := a\a k . 
The Hamiltonian of a free phonon field reads 

(2.4) tf ph = I dkuj{k)N k . 
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Here uj(k) is the dispersion relation for phonons. The state of thermal equilibrium for a phonon field 
at inverse temperature (3 and chemical potential p is given by the density operator 

(2.5) 7ph :=Z- X e -l3H ph +^jN k dk ; with z ._ Tin ^^ e -pH ph+f ifN k dky 

The expected number of phonons in the mode k reads 

„— /3u>(k)+n 
M(k):=Tr nph ( lph N k ) = — — 

The symbol C will denote various positive constants arising in estimates. 

3. Assumptions 

In the following we will use extensively the notation (x) := (x 2 + l)a and the estimates 

(x + y) <C(x)(y) 
(3 ' 1} (x + y)" 1 ^Cixy'iy). 

The dispersion relations e(-), cj(-) for electrons and phonons are assumed to be: 

• spherically symmetric; 

• decaying for large k up to the 2cP h derivative: 

|V J e(fc)| < C(l + (k) 2 - 1 ), I = 0,1,2,..., 2d; 
\V l u(k)\ < C(l + (k) 2 - 1 ), I = 0,1, 2,..., 2d. 
We will often consider the functions 

(3.3) $±(p,fc) := e(k+p) ±u{k). 
We assume that the following relations hold 

(3.4) lim $±(p,k) = oo 

fe^oo 

(3.5) 3 Ci, C 2 £l s.t. < d < Hess fc $±(p, k) < C 2 . 
Let us consider the "thick level sets" 

E±(p,d,5) := {k : \*±(p,k) - 6\ < 6} 

for the functions <&±. From (|3.2p . (|3.4p . (|3.5p one can conclude that there exist C, p such that for any 
8,p < p 

(3.6) sup\E±{p,e,5)nB{q,p)\ <C5p d - 1 , 

where | • | denotes the Lebesgue measure, and B(q,p) is the ball of radius p centered in q e We 
require a further condition on the intersection of two such sets: there exist p, C3 > such that for any 
5i,5 2 ,p<p, Pu p 2 £R d , 0i,02 GK 

/T x x d 2 

(3.7) sup\E±(p 1 ,9 1 ,6 1 )nE ± (p 2 ,e 2 ,6< 2 )nB(q,p)\ < , 3 1 2p . 

q \Pl-P2\ 

The listed conditions are fulfilled for a classical or a pseudo-relativistic electron if HV^Hoo is sufficiently 
small. 
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To ensure that the density operator for the phonon field is trace class we assume 

inf w(Jfe) - M/?" 1 >C>0. 

k 

The electron and the phonon field will be coupled via an interaction form factor F(k) : M. d 
assumed to be real and symmetric, namely 



It is 



F(k) = F(-k) = F(k) 
and it has a fast decay up to 2d-derivatives 



(3.8) 



max 
=0,...,2<2 



V{F(k) < C{k) 



-2d- 12 



For a = ±1 we define the functions 



L(k,a) := |F(fc)| 2 (M(k) + 



a + 1 



and 
(3.9) 



£{k) := L{k, 1) = \F(k)\ 2 (M{k) + 1) . 



Obviously, L(k, ±1) < C(k), and C shares with F the decay estimate f|3 . 8|) . 

These conditions on e(k), u(k) and F(k) ensure the following key estimates that we will subsequently 
use: 

Lemma 3.1. For the functions C defined in t3.9\) , and &± defined in \3. 3\) , the following estimates 
hold 



(3.10) 
(3.11) 
(3.12) 



sup / 

p,e Jm 



C(k) dk 



J 



9-$±(p,k) +in\' 
C{k) dk 



- < Cn' m , m > 0; 
C(log^) 2 



\9 - $±(f>, k) + in\ \8 - $±(w, k) - irj\(p)(u) 

da 

sup 



< 



\p - uU(9)2 (9)2 
< Clog* 7?; 



\a - e(v 2 + i/2) + irj\ (a) 
where rj > and we used the notation f* = max(l, |/|), f* = min(l, |/| +77). 



The new estimate compared with [E] is (|3.1ip . whose proof is given in the appendix. 
In the proof of Theorem 11.11 we will use the following inequality too: 



(3.13) 



sup 

V 



E 

a G {±} 



< 



c 



The proof of (|3.13p is in Lemma 4.1 of [Ej. 
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4. The model 

4.1. Hamiltonian. The dynamics of the system consisting in an electron interacting with a phonon 
field is generated by the Hamiltonian 

H = H e ®I nph +I Hc ®# ph + tf/. 

In any tensor product the first factor acts on TC e and the second on TC p h- The free Hamiltonians H e 
and -ffph have been defined in (12. ip and (12. 4|) . The interaction Hamiltonian Hj is defined by 



(4.1) Hi = iXJ dkF(k)(e- tk - x al-e lk - x a k ). 

As stated in Section [1] we will always assume A = yfi. 

We assume that initially the electron and the phonons are independent and their initial state is 
represented by the density operator 

To := 7e,o ® 7 P h 

where 7e,o is an electronic density matrix and 7 p h is the phonon thermal state defined in (|2,5p . The 
time evolution of T t is described by 

idtT t = [H,T t ] 

with Tq as initial data, or, equivalently 

1 1 — e 1 oe 

4.2. Initial state of the electron. The initial state of the electron is represented by the wave 
function (jl.ip . The density operator corresponding to the wave function i/jq is the one-dimensional 
projection 

7e,0 = |-0g> <-0§ I 

with integral kernel 



According to (|l.ip . the initial density operator for the electron can be split in four terms 

( 4 - 2 ) 7e,0 = 7e,0,++ + 7e,Q,— + 7e,0,+- + 7e,Q,-+ 

where 

(4-3) 7e,0,aa' = iV'o.a) ($),<*' I 

with a, a' £ {+, — }. We stress once again that, whereas 7 ej o,++ and 7^0, describe an electron in the 

state ipQ + and ipQ_, respectively, the terms 7 e ,o,H and 7 e ,o, h do not represent any physical state: 

in particular, as operators they are not positive. They represent the interference between the states 
V>o,+ and 

According to (14. 2p , the Wigner transform of the initial state for the electron can be written as 



(4.4) W % = W ,++ + Wo,— + W ,+- + W 



RATE OF DECOHERENCE FOR AN ELECTRON WEAKLY COUPLED TO A PHONON GAS 



where 



W , ++ (x,v) = W f {ex + Q, V —^j 

W 0t —(x,v)= W f (sx-Q,^^j 

W 0l+ -(x,v) = e 2lP - x ecW f (ex,~ 

W 0t . + (x,v)= e- 2ip - x e- 2j ^W f (Ex, V - 



as one can immediately verify by (|2.2p . 
Recalling the convention adopted in (|2.3 



<4 ' 5) f ^.. ) . ( v &) W(= + ^)?(=-4r)- 

All components of Wo depend on e, but this is omitted in notation. 

4.3. Time evolution of the electron. The state of the electron at time t can be described by the 
marginal (or reduced density operator) 

(4-6) 7 e,t := Tr Hph r t 

or, equivalently, by the Wigner function 



Since the time evolution is linear, the initial decompositions (14. 2ft and (14. 4h propagate at any time t: 

7e,t = 7e,i,++ + 7e,t, + 7e,t,+- + 7e,t ,-+ 

(4 ' 7) W(t) = W++(t) + W—(t) + W+_(i) + W-+(t) 

Again, the first two terms in the decomposition are called diagonal, and the other two terms non- 
diagonal. 

4.4. Observables. We will use the notation O for the observables considered as operators and their 
Wigner transforms will be denoted by J = Wo- The expectation value of an observable is given by 



Tr C*7c = / J(x,v)W(x,v)dxdv = I J £ (u,v)W (u,v) dudv. 

We need observables capable to resolve the size of the interference fringes. Therefore we consider 
observables O = O e endowed with a two-scale structure in the space variables but we always assume 
the condition ([1.5)1 on J £ . This implies that 

(4.8) |||0||| 2 := sup||0*0 e || < oo 

e 

since the operator norm of 0*0 is bounded by || J|| 2 . 
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In order to select observables that actually detect the exponential decay of the interference fringes, 

we provide sufficient conditions on J £ for lim e _»o(^e, i f tee (e~ 1 T)) to exist. We assume that the 

Wigner transform of O e has the structure 

(4.9) Je(x,v) = A(ex,v)b(x) 

where A describes the macroscopic profile of the observable and b encodes the short-scale structure. 

l 

2 

V 



First we consider for simplicity the case of Newtonian dispersion relation, i.e. e(v) = \v 2 . One has 



W+^fU^T; x, v) = e 2ip -*e 2i T^- pT )W f ex - vT, - 



£ 



Due to the regularity hypothesis ([T2]) . W f belongs to L 1 ^ 2 ^), and then ||W+_ f ree (e -1 T)||i = ||W/||i. 
Besides, due to (|1.5p . it is clear that J £ £ L°°(M. 2d ) with a uniform bound in e. 
After an elementary change of variables, one gets 

(«/ £ ,^+-,frce(^ 1 T)) = j dxdvA{x + evT,ev)b g + wT) e 2iP i^ +vT )e 2lv < Q - pT ^W f (x,v) 

and the limit for e — > can be investigated using the dominated convergence theorem, assuming that 
A and b are bounded and continuous functions. One finds that: 

• if liniM^oQ b(x) = 0, then lim e ^o(^ei ,free( e ~ 1 ^ 1 )) vanishes. In this case the observable 

does not resolve the interference fringes, even though the limit of the expectation value exists; 

• if lim^i^oo b{x)e~ 2lP x =: q, ^ 0, then 

(4.10) lim(J £ ,W + _ ifree (e- 1 T)) = ~ h J dx dv A(x, 0)e 2iv < Q - pT) W f (x, v) 

In general, the short-scale factor b must exhibit the periodicity of the interference fringes. We 
require the distributional Fourier transform &(£) to be a complex measure with a finite total variation. 
If it has a non-trivial Dirac delta component with frequency 2P, i.e. 

(4.H) lim [ b(0dti=:c b ^Q, 

then the interference fringes in can be detected by J e . For the other off-diagonal term, W |_, 

one needs a non-trivial component with frequency — 2P. Working in Fourier space requires more 
conditions on A, for simplicity we assume that A is a Schwartz function. 

In the case of a more general dispersion relation e(v) the Wigner transform of the free evolution, 
W+-,free(*), is given by 



(4.12) W+^mv) = (^ e )- d e ^^/ 2 )-e(,- 5 / 2 )] e 2».f + £Z|^ f(l - LJ^ 

in Fourier space. After a change of variable one gets 

(J e , W+^beeie^T)) = J d£ dv d(I U + Sv) %) e -^-'[e( CT+ ^/2 + P)- e ( CT - £ e/2-P)] 



xe 2i ^f{v + i)f(v-i 
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Again, by using dominated convergence, we have 



ton(J e ,W + -j M {e- 1 T)) = (2ir)- d /% J d£dv %0) e 2 ^~ TVe ^))/ (v + f(v - |) 

that, recalling (|2.3p . generalizes (|4.10p to a generic dispersion relation. 
In summarizing, we have proved the following 

Corollary 4.1. Assume the hypotheses of Theorem \l.l\ and let J e (x,v) = A(ex,v)b(x). Let either 

i) e(v) be quadratic, A £ C b (R 2d ), b £ C b (R d ) and c b := \\m\ x \^ 00 b(x)e^ 2iP ' x ; or 

ii) A € 5(R M ) and b be a measure with finite total variation with c b being the coefficient of the 
Dirac delta component at 2P of this measure (see (I4.11j) ). 

Then, 

lim lim<J e ,W-+_(e- 1 T)) = e~ T<Jp lim( J e , W^^e^T)) 

= e~ T ^- b J dxdvA(x~^)e 2m < Q - TVe ^W f (2 

The hypotheses cover the case J £ (x,v) = e 2iPx for the quadratic dispersion relation (see (|1.3j) ). 

5. Known results 

The proof partially follows the arguments in [E). To avoid duplicating them but still keep the 
present paper relatively self-contained, here we summarize some notations and results from [E] . 



[x,v 



5.1. The main term of 7 ej t. We consider the density matrix of the whole system 

Lt — e 1 oe 

and expand the propagator e~ lHt into a Duhamel sum up to the order 



log I loge| 

After performing the phonon trace and defining the reduced density matrix of the electron as 

7e,t = lr ph (e 1 oe ) , 
7 ej t is decomposed into the sum of a main term and an error term 

(5-1) 7e,t = 7r n (0+7FW 

where K £ N is an arbitrary integer. For the precise definitions of 7j£ am (£) and 7^ r (t) see formula 
(2.18) in [E], but their explicit form will not be needed here. We need only to know that 

(5.2) lim sup lim sup lim sup | Tr 7f r (e _1 T)0*| = 0, 

K— >oo e— >0 L^oo 

i.e. that the estimate (2.24) of [EJ holds in our case as well. To arrive at this result, the argument 
in [E] used only assumption (1.26) in [E] on the initial state and assumption (2.21) in [E] on the 
observable. In our case, the first assumption is guaranteed by the hypothesis fjl.2[) ; the second one 
is inequality (|4.8p . For the same reason, Lemma 3.1, Proposition 10.1 and Lemma 10.2 from [EJ also 
hold in our case. 
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To evaluate 7]^ am (i) we use Proposition 10.1 and Lemma 10.2 in [EJ. Results therein can be resumed 
as follows 



lim sup 

£^0 



(5.3) 



K-l min(AT,jV) 

E E 



N,N=0 n=0 m = (mo, . . . ,m„) e N" +1 m = (m , . . . , rh n ) e N n+1 

where, for any m, m G N™ +1 , i > we defined 



m,m, id Vv 



^e{±} ^ \i=l / V i=l / 

j=l,...,n 



d5e it5 m^ +1 T^ 



(5.4) 



Here 77 > 0, k = (ki, . . . , k n ) G W 1 , |m| = X]j=o m^, |m| = X^™ =0 fhj and to define i?j and Tj we first 
introduce 



R(a, v, z) 



1 



(5.5) 



a — e(v) — z 



T v (a,v) = / - 

<T = ± 17 



L(fc, cr)dA; 



e(u + k) — cru(k) + irj 



for a £ I, D £ R , z G C, and 77 > 0, and then we set 



R i 



Rla,v+ ^2 fc i + 9> X] ^M^) 



1=3+1 1=3+1 
n n 



Tj := T,, a- ^ aiuj(k),v+ ^ fej 



RATE OF DECOHERENCE FOR AN ELECTRON 



WEAKLY COUPLED TO A PHONON GAS 



13 



Formula (|5.4j) holds for any rj > 0; later we will choose r] = e. 
The following estimates will be used in the next section 

sup IT^ (a,p)\ < C, \Rj\ m i +1 < v~ mj \Rj\, 

(5 6) P ' a,V 

sup (\V p T v (a,p)\ + \d a T v (a,p)\ + |$,T„(a,p)|) < C^ 2 

The first and the third inequalities in (|5.6p are proven in [E] (Lemma 4.1), the second is trivial. 

5.2. Ladder Wigner functions and related Feynman graphs. We introduce the family of "lad- 
der Wigner functions" W^^^t), defined implicitly by the relation 

The explicit expression of W^ d ~ (t) in the Fourier space can be easily obtained comparing (|5.7p with 
(|5.4p . The content of (|5.3p can be expressed saying that, in the thermodynamic limit, the expectation 
value of O s on 7^ am (t) can be decomposed in the sum of the expectation values of J £ on the states 
represented by W^g(t). 

The function W^^f^fa^v) is linear in W^, thus it contains two evolved copies of the electron 

wavefunction: one copy initially lies in the state represented by ipQ, the other in ip^. In the following 
we will call them, for short, the first and the second copy. In the time interval (0, t) each copy of the 
electron emits and absorbs phonons. An important feature of the evolution described by W^^^t; £, v ) 
is that any emitted phonon must be reabsorbed and vice versa, but a phonon emitted by one copy of 
the electron can be reabsorbed by the other. 

This fact permits us to distinguish two different processes: the exchange of phonons between the 
two copies of the electron, and the recollision of the same copy with the same phonon. In the evolution 
described by W^^^t; £, v) all possible recollisions are immediate, i.e., there are no further events in 
the time interval between the emission of a phonon by a copy of the electron and its reabsorption by 
the same copy. Moreover, exchanged phonons appear in the same order in the time evolution of the 
two copies of the wavefunction of the electron. 

In (|5.4|) there are n exchanged phonons and their momenta are labeled k\, . . . ,k n . The electron 
propagator between consecutive phonon collisions is represented by the factors Rj (for the first copy) 
and Rj (for the second). Between the jth and the (j + l)st exchanges of phonons, the first copy 
of the electron undergoes rrij immediate recollisions, each of them embodied in a factor Yj in (|5.4p . 
Analogously, the second copy of the electron suffers fhj immediate recollisions, each of them represented 
in (|5.4p by a factor Yj. 

The expectation value of an observable O e on the state represented by W^^^fa^v) can be rep- 
resented as a Feynman graph (see Fig. 1). 

5.3. Non-diagonal terms. Following (|4.7|) . the expectation of any ladder Wigner function can be 
decomposed in diagonal and non diagonal terms 

^m,rh,id (^) ^ ] ^ aa! , m , m ,id (^) 

a,a'e{+,-} 
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k n ~2 
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k n 


n n 




n n n 




n 



«-e/2 



Figure 1 . A ladder Feynman graph. The two horizontal lines represent the time evolution of 
the two copies of the electron, and the vertical lines joining the two horizontal ones represent 
the exchanges of phonons between them. They are labeled by the momenta of the exchanged 
phonons. By convention time goes from the right to the left and the labeling of the exchanged 
phonons follows the opposite ordering. The square loops on the higher horizontal line sym- 
bolize the recollisions in the evolution of the first copy of the electron. Such loops have no 
momentum labeling because in the related factors Tj an internal loop momentum integral has 
been performed. An analogous interpretation holds for the square loops in the lower horizontal 
line. Here m = 1, m„_2 = 0, ?n n -i = 5, m n = 2 and m = 2, m„_2 = 3, rh n -\ = 0, fh n = 1. 



1>< 



where C* a , m mid(*) * s defined from (|5.4p after replacing W^g by Wo^a'- The corresponding Wigner 
transforms, W^£,^~(t), are denned implicitly via 



(5.8) 

and of course, 



C-\ — ,m,m,id(^) — (^S) ^aa' ,n,m,m(^)) 
gladder = 



ETj/ladder t±\ 
act' ,n, m , m \ ) 

a,a'e{+,-} 



Since both the evolution of density matrices and the definition of 7j£ am (t) and 7fp"(i) are linear 
in the initial data (see (2.18) and (2.19) in [E]), all results listed in Section 15.11 hold separately for 
diagonal and non-diagonal terms, according to decompositions (14. 7p . In particular, formulas (|5.1j) . 
(|5.2p . (|5.3p remain valid after replacing 7^, 7^ am (i) and 7fp(i) with the corresponding non-diagonal 
component 7 e ,t,H , 7^ a | n _(t) and 7fJ r H (t), respectively. 

Then, formulas (|5.2p . (15.3H . applied to non-diagonal terms give 



lim sup lim sup 



lim (J £ ,W+ 

L— >oo 



.-1 



T)> 



K—l min(JV,JV) 

E E 



E 



E 



C* 



,m,id( e T) 



N,N=0 n=0 m = (mo,...,m„)eN" +1 m = (mo, . . . , m n ) G N n+1 
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In [E] it is proved that the limit 



lim lim lim [TV 7^ ain (e~ 1 T)C'* 

K^oo e^O L—*oo 



equals the expectation value of O e on a distribution on the phase space that solves a linear Boltzmann 
equation, provided that O e has no short-scale structure (i.e. b = 1 in (|4.9p ). 

In the rest of the paper we study the same limit of the non-diagonal part for a two-scale observable 
(|4.9p . In the following section we prove that ladder terms with at least one rung do not contribute. In 
contrast, diagrams with no rungs and pure recollisions give a non vanishing contribution, characterized 
by the exponential damping described in Section [TJ This is the content of Section [7J 



6. Non-diagonal ?i-rung ladder term with immediate recollisions: estimate 

We estimate the contribution of the generic non-diagonal ladder term W x ^ A ^ m ~ (i) to the expec- 
tation value of a two-scale observable J £ . 

Lemma 6.1. Given n > 1, m,m e N n+1 , P / 0, let W^ dd ^ mfh be defined as in (|5.8fl . and J £ be a 
two-scale observable satisfying estimate (jl.5p . Then for any fixed T > 



(6.1) (Je,W^ d * Ae- l T)) < C"+H+^||J||^(log* £ ) 4 £ 



\P\ 



H 2 (K. d )i 



where \\J\\ is defined in (|1.5p . |P|* = min(l, \P\ + e), log*e = max(l, | log e\) and the constant C 
depends on T. 

Proof. From (|43j) . ([OP and (|51?P we have 

(J £ , W_\ — : n,m,jrSS)) 

= (2tt)- 2 -1 e- d \ 2n+2 ^ +2 ^\ Y, [d(dv (f[dkA J £ L,v + j2k)e 2lV -f 

CTj e{±} J \j=l J \ 1=1 J 

j=l,...,n 

Vj=i 
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To estimate \{J S , W^_^ m ^(e T)) | we use the first two inequalities in (|5.6|) . with rj = e, and definition 
(GDI. Then 



' 1 £ 2e 



• 1 e 2e 



j=l,...,n 



(6.2) 



Let us focus on the integral in dk\ . . . dk n dada. 



f{J C{kj)dk}j Jja (j[\Rj\J 



da Y\ I Rj I 

3=1 



(6.3) 



dada\R n \\R n \ / dk n C(k n )\R n -i\\R n -i\ . . . \ dk\ H{k\)\R^R 



Indeed, the variable k\ appears in Rj and Rj only if j < I — 1. Therefore, we can integrate in (|6.3p 
following the order k±, . . . , k n . The integrals in fc,- with 1 < j < n — 1 are estimated by standard 
Cauchy-Schwarz inequality 



J Cikj^Rj-iWRj^ldkj < C J £(^)l^i-i| 2 dkj + C J C{kj)\Rj^\ 2 dkj < Ce~ l 



(6.4) 



where we exploited estimate (|3.10j) . with rj = e and rn > 1. We finally obtain the factor C n 1 e n+1 . 

The integral in k n is estimated using (|3. with 6 = a, 9 = 5, p = v + £/2, p = u — £/2, rj = e. We 
obtain 



£,(k n ) dk ri 



\a - e(v + f /2 + fc„ 



a n uj(k n ) + ie\\a - e(v - £/2 + fc„) - a n uj(k n ) - ie\ 

f (log*e) 2 <T, + £/2><i;-£/2> 



< C- 



|eU(a)V2(5)i/2 



(6.5) 



where, according to the notation used in Lemma [37TI we defined = min(l, |£| + e). 
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Finally, by (I3.12p both integrals in da and da give a factor log* e. Then, by (|6,2p . (16.41) and (16.5 

we have 



< e- a {CX z e 



2 ~l\n+\m\ + \m\ 



(log*e) 4 e / d£dv 



' 1 e 2e 



f 



v j-2P 
e 2e~ 



leu 



(6.6) 

We focus on the integral factor. Using (|1.5p we get 

d£dv J e (f, 



(6.7) 



< ||J||sup 



1 



1 e 2e 



1 e 2e 
1 e 2e 



{v + i/2){v-i/2) 

leu 

( V +e/2)(«-e/2) 



#2( K d) SUp 



1 



sup 



<t; + e/2)<w-C/2> 



/r , g-2P \ 2 /i, g-2P \ 2 
2e / \s 2e / 



+ 



By (|3.ip one has 
(6.8) 

and 
(6.9) 

Elementary calculations show 
(6.10) 



( V + £/2>(«-£/2) 



v , g-2P 
e "•" 2e 



e 2e 



« g-2P \ Iv j-2P \ > 
e 2e / \ e 2e / ~ 



£-2P \ 
2e / 



sup , 



< 



c 



T5-/iei 

Then, plugging (|5^]> . (IBT91) and (IBTTUD into ((SZZD we obtain 

'v £-2P 



<i£ dv 



1 e 2e 



<c e W P) ' 



2s 



{v + Z/2){v-Z/2) 

leu 



LP 



□ 



By (|6.6p . recalling that A = we conclude the proof. 

Remark 6.2. For n > 2 one can gain another e factor in estimate f)6. 1 1) . This is easily accomplished 
using inequality 

L{k) dk C(log* r]) 2 



sup 



< 



|0 - $±(p, fc) + ^| |0 - *±(«, fe) - "?l b - n l 

for estimating the integral in k\ in (16. 3p . 
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v + C/2 - v + £/2 

□ □ □ □ □ □ □ " 



4' 



e 



□ □ □ □ $ 

Figure 2. An example of Feynman graph associated to a pure recollision term. Here m = 7, 
m — 4. 



e 



7. Terms with immediate recollisions only 

In this section we prove Theorem 1.1, i.e. we compute the contribution of terms consisting of 
immediate recollisions only, namely terms with n = in f)5.3[) . The corresponding Feynman diagrams 
are illustrated in Fig. 2. We stress that, according to the results in the previous sections, these are 
the sole non vanishing terms. 



Proof. Since the vectors m and m are one-dimensional, we can simply denote by m both m and its 
only component uiq, and by fh both m and mo- 

We preliminarily observe that, according to (|5.4p . (|5.7p . (|4.5p . and the definition of non-diagonal 
terms given in Section [5T3l for any t > and any r\ > the function Wrr^ rofi (f) can be expressed 
as follows: 



V £" 2S" y V <S 2£" y 

x ( dae- ita R m+1 (a,v + £/2,-i7])r™(a,v + I da e itS i?™ +1 (5, v - f/2, irj) 

Jr ' Jr 

xT™(a,v-{;/2) 



(7.1) 
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A 



2m+2m 



-d 



2vr £- a e 



d„2i- 



v_q 7 (v f-2P\ 7 (v i-2P 



f . (2m 

/ n< 



7 u + 

2m 



2s 



2s 



Hd Sj \6 h-J> [U e 



■is 2j e(v+(,/2) 



j=0 



n e 

j=i \<re{±} 



dfe L(k, a) e -^-i[e(u+£/2+fc)+™(fc)] 



2m 



is 2 je(v-(,/2) 



j=0 



U"=0 



n e 

J = l \<re{±} 



-i»2j_l[e(n-^/2+fe)+crw(fc)] 



(7.2) 

where the functions R and were defined in (|5.5p . The expression at the r.h.s of (|7,ip is said the 
resolvent form, the one in (|7.2p is said the propagator form of the ladder Wigner function. One can 
easily obtain the first from the second one using the identity 



2m 
3=0 



2^ 



-iat e i(a+ir,)Y: 2 ^o^ da, 



integrating in the variables Sj, and proceeding analogously for the variables Sj. 
Moreover, we introduce the following quantity: 

L(k, a) dk 



(7.3) 



:= lim > / 

v-^o+ J 



e(P) - e(P + k) — aoj(k) + ir] 



T 0+ (e(P),P). 



The existence of the limit has been proven in Lemma 4.1 of [Ej. 

Besides, from the definition of the Boltzmann collision kernel (see (1.27) in [E]), the total cross 
section ap for an electron with momentum P in the phonon field reads 

(7.4) a P = 2ir ^ j dU L(P - U)5(e(P) - e(U) - au;(P - U)) 

and identity (4.12) in [E] yields 

Im <3?p = 

2 

We organize the proof into four steps. 



(7.5) 



Step 1. Here we prove 
(7.6) 



< 



{C\H) 



2i\m+m 



ml ml 



Let us mention that this estimate is stronger than (10.4) in [Ej, but applies only to the ladder terms. 
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We use the exponential representation of the propagator (I7.2p . Integrating the delta function with 
respect to the variable S2m one gets 



J[0,t]2m+l \JJ- J 



2m 



6 t-^sA R, 



fl E / dkL(k,a)e 



3=0 / V=0 

is 2 j-i [e{v+k+(,/2)+aw{k)\ 



j=i We{±} 



ft ft — S() 

Jo Jo 



dso. . . 



*-£ i=0 s 2j 



2m-2 



*-E,-0 



dsi ( ^ / dk L(k,a)e- ls ^ v+ V 2+ V + ™W-< v+ tm] 



d«2m-i [ 2 dkL(k,a)e- tS2 ^< v+ ^ 2+k ^^-< v+ ^ 



Now we can use inequality (|3.13p and obtain 



3=0 J \j=0 



-is 2j e(v+£/2) 



f (fids a 

J[o,t] 2m+1 \fJo J 

f[( Y, [dkL(h,a)e- is v-M v+ V 2+ Q + ™W J 
j=i Ve{±} J ) 



< \ ds 

'0 JQ 



t-s 



dso... 



*-Ej=0 s 2j 



2m-2 



U = l 



n / ds2 J- 1 



Since the space dimension is at least three, the quantity at the r.h.s. can be estimated by 



C m / ds 
'o JO 



t-so 



ds 



2 • 



*-Z J j=0 s 2j C m t m 
dS2m-2 = — 







An analogous estimate holds for the factor in (|7.2p involving the variables sq, . . . , s„. 

v i-2P 



(Je, ^+ a - d o e ,m„m(*)) 



{ i- e~ d I d£dv\J E (Z,v) 



m\m\ 



< 



(CXH) 



2+\m+m 



m\m\ 



(f£ sup 



£-2P 



/ 



dv 



f v + 



2s 



e 2e 

i-2P 



f v 



2s 



where we performed the change of variables v — ► sv. Estimating the integral in v by Cauchy-Schwarz 
inequality and using (jl.5p we prove (|7.6h . 
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Step 2. Here we prove 



(7.7) 
where 



lim 



ladder 



-2—d/2 —d \ 2m+2m ff.ru jwmjli - 



v t £ - 2P\ j(v i-2P 



2s 



■ 



2s 



:e 2tri [ dae- ita R m+1 (a,v + Z/2,-iri) [ da e ita R m+1 {a, v - £/2, irf) 



(7.8) 

and $p was denned in (|7.3j) . Comparing formula (|7.8p with (17. ip . one easily notices that W 1 ^^ ~ 
equals W^^^^ apart from the replacement of by $p. Then, in the case m = fh = there is 
nothing to prove. We give the proof in detail for the generic case with both m and fh greater than 
zero, and leave to the reader the case in which either m or fh vanishes. We have 

= ^y 2 ^/ 2 s- d \ 2m+2ih e 2 ^ j didvlA^v)- 2i — 
< / dadae- ita e i& R m+l R fh+l 



' 1 s 2s 



v _ £- 2P\ j(v £-2P 



2s 



(7.9) 



where for brevity we used the notations 

R = R(a,v + £/2,-i V ), T v = T v (a,v + £/2) 

R = R(5,v-£/2,iri), T v = T„(5, v - £/2) . 
By the first estimate in (|5.6p we have 

fh— 1 m— 1 



T- (T, - Ip) £ Tj S^'" 1 + *5 (T„ - * P ) £ Tj cDp^- 1 



j=0 



(7.10) < C7 m+ ™(|T r/ -$p| + |f^-¥p 
Now we estimate (|7,10p . From the third formula in (|5.6p one has 

|T„-$ P | < Ci 1 - 1 ' 2 {\a-e{P)\ + \v + t/2-P\+r ] ) 

(7.11) < Ci 1 - 1 / 2 {\a-e{v + i/2)\ + \e{v + i/2)-e{P)\ + \v + i/2-P\+r 1 ) 
Moreover, the first estimate in (|5.6|) shows that 

(7.12) \r v -<f> P \ < \r v \ + \$ P \ < c 

From (fTTTj) and (T?Tl2|) we obtain 

(7.13) |T„-$ P | < Cr ? - 1 / 2 (|e(v + ^/2) - e(P)\ + \ v + £/2-P\ + min(|a - e(u + £/2)|, 1) ) 



22 



RICCARDO ADAMI AND LASZLO ERDOS 



where we used rj = e < 1. In order to obtain an analogous estimate for 
that, due to the symmetry of the functions e(-), u>(-), L(-,cr), one has 

¥ P = T 0+ (e(P), P) = T 0+ (e(P), -P). 

Hence, mimicking (|7.1ip . (|7.12p . (|7.13p . one obtains 



we first observe 



(7.14) 



r v -<s> P 



< Cr]~ l / 2 {\e{v - f/2) - e(P)| + \v- f/2 + P\+ min (\a - e(v - £/2)|, 1) ) 



From (JE2I), (HHU]), (|715|> . ([731D we obtain 

(J^W^^ie-'T)) - (J £ ,W^ m ^(e- l T)) 



d —1/2 //^< \ 2\m+m 



d£ (it" 



/ 



1 e 2e 



u £ - 2P 

£ 



2e 



da da \R 



m+l I p|m+l 



\R\' 



\e{v + t/2)-e{P)\ + \e{v-i/2)-e{P)\ 



+\v + i/2-P\ + \v-i/2-P\ + min(|a-e(i; + £/2)|,l) + min (|5 - e(v - £/2)|, 1) 
(/) + (/') + (//) + (//') + (///) + (///') 



(7.15) 



where the last decomposition is made according to the terms in square parentheses in the r.h.s. of 
(171511 . 

The first term is estimated as follows 



(I) < e- d r)- 1/2 {C\ 2 ) m+ih J dCm V \J £ (^v)\ I dv 

x|e(u + f/2)-e(P)| f dada\R\ m+1 \Rf +1 
Recalling that m > 1, we have 

da: 



1 e 2e 



f 



v £-2P 
e 2e~ 



\a-e{v + i/2) +ir]\ m+1 



< Cif 



and analogously for the integral in a. 

By Cauchy-Schwarz the integral in v gives 



dv 



1 e 2e 



1 e 2e 



| e ( u + £/2)-e(P)| 



<i/2 



1 e 2e 



| e (t, + £/2)-e(P)| J 



Changing variable to u := 7 + ^ — 7" ° ne gains a further factor e d//2 . Besides, using hypothesis 



3.2p one gets 



\e(v + t/2)-e(P) 



\e{£u + P) - e{P)\ < C£(u) 2 {P) 



RATE OF DECOHERENCE FOR AN ELECTRON 



WEAKLY COUPLED TO A PHONON GAS 



23 



Then 



dv 



1 e 2s 



f 



v £ - 2P 



\e(v + i/2)-e(P)\ < Ce d+1 (P)\\f\\ 



e 2e 

After getting rid of the integral in the variable £ by (jl.5f) we obtain 

(/) < er lf \c\*r l ) m+ih {p)\\f\\HH^) < c^sy 2 (p)\\f\\ H , 



where we chose rj = e and used A = \J~k~. 

Term (II) in (|7.15p has the same structure as (/), apart from the replacement of e(v + £/2) — e(P) 
with v + i/2- P. 

For estimating (III) we proceed as in lemma 10.3 in [Ej. Using the notation 9(s) = min(|s|, 1), we 
observe that 



(a-e(v + t/2)) 



< 



C 



\a - e(v + f/2) + iri\ ~ {a - e(v + f/2)) 

Then, by assumption (11.51) . hypothesis m > 1, Cauchy-Schwarz inequality and the following two 
estimates, easily derived from (|3.10p . 

da 



{a - e(v + f/2))|a - e(v + £/2) + mj|» 

da 



|5-e(u-£/2) + i7/| m+1 



we conclude 



(III) < e- d (CX 2 ) m+m r]- 1/2 J d^dv \j(£,v) 



/(^" 2P 



e 2e 



1 e 2e 



da da 



(a - e(u + £/2))|a - e (v + £/2) +ir/| m |a - e(v - £/2) + zr/| m + 1 

For terms (/'), (II 1 ), (HI 1 ) in (|7.15p we have the same estimates found for the corresponding terms 
without prime. Then, 

(Je,W^ m ^(e' l T)) - (Js,W l ^ m ^(E- l T))\ < CT**y/i(F)\\f\\ m 

and (|7.7p is proven. 
Step 3. Here we prove 



(7.16) (Je,y^% n , iA {e- 1 T)) 



-iT$ P ) m (iT$ P ) 



ml 



ml 



-(J E ,W + < Tee (e- l T)) 



First, notice that the integrals in da and da in the definiton (|7.8p of W^f_Q r m ^(t) can be explicitly 
computed using residue theorem. One obtains 



-ita 



da 



[a - e(v + £/2) +irj] m+1 



-2me~ m e 



nt-ite(v+t/2) 



(-ay 



m\ 
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and analogously for the integral in da. Therefore, choosing 7] = e, t = e 1 T, A = ^/i, from (|7.8p one 
obtains 



{Je,yv+-« m ,m[£ T)) = £ (2vr) 



m 



m 



d£dv e^TW^/ 2 )-^-?/ 2 )] 



xj e (£,v)e 2t e f _ + s * < ' 



2s 



2s 



(7.17) 

To compute JV H free( e_ T) we write 



7e )E -iT 1+ -,free = ^1^(11^ 



or equivalently, in the Fourier space 

7e, £ -lT,+~,free(P)«) = 



| [e(p)-e(«)] e i(p+«)-f f^l£)j j 



U + P 



Finally, (|7.16p immediately follows from (17.17P by using (14.121) . 

Step 4- Here we prove (|1.6p . By formulas (|5.2p . (|5.3p . (|5.4p . (|5.7p for non-diagonal terms, and 
eliminating terms with n > by Lemma 16-H one gets 



lim 



lim(J e ,W f _(e- 1 r)) - e- T ^(J e ,Ty + _, frce (e- 1 T)) 



L—>oo 



lim lim 



m,m=0 

where we used (17. 6h to exchange the limits in K and in e. Then, from (|7.7p and (|7.16p 



lim 

£^0 



lim (J e ,W + _(e- 1 T)) - e- r ^(J e ,W + _ ifree ( £ - 1 T)) 



L— >oo 



lim lim 



< lim sup 

K— >oo 



m,m=0 



-Tap 



ml 



m! 



ml 



ml 



limsup (J e ,W+_free(e T) 



m,m=0 

The first factor in the r.h.s. vanishes due to (|7.5p . while the second factor is bounded since 



\(J e , W + ^ hee (t))\ = Tr O^e-^-^+Xe-**^^. 
The proof is complete. 



< \\o. 



□ 
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8. Appendix 



Here we prove Lemma 13. 11 

Proof. Inequality (|3.10p corresponds to formula (5.19) in [EJ, and we refer to that paper for the proof. 
We prove (|3.1ip . We treat the case of since for the proof is the same. 

Consider a tiling of the A:-space in cubes {Qi}igt*j of size thinner than p defined in (|3,6p . (|3.7p . Fix 
Iq := [log* (pr/ -1 )], and define the sets 



S, 



~-{k e 
--{k e 
~-{k £ 



e~ l ~ l p <\9- $+(p, fe)| < e' l p}, 1<1<1 -1 

\9-<5> + (p,k)\<e- lo p} 

\9-$ + (p,k)\>p/e} 



and analogously the sets Sj, < I < Iq, where 9 replaces 9 in the definitions above. 
Then, 



C{k)dk 



lo 



\9 - <f>+(p, k) + irj\ \9 - <3?+(u, k) - irj\(p}(u) 



1 1,1=0 



where by | • | we denoted the Lebesgue measure in M. d and introduced the notation 

( C(k) \ 



U. ,r:= 



sup 



keQinS^Sy \\0- $+(p, k) + irj\ \9 - k) - vq\ (p)(u) 



Since the actual value of p does not play any role, in what follows we will absorb it in the constant C. 
We prove the following estimate: 

(8.2) U ilT < C sup[(fc) 4 |£(fc)|], 0<l,l<l o . 

{9) 2(B) 2 k 

First, we treat the case 1 < 1,1 < Iq. We have 



U. , r < sup - 

k£ Qt n Sl ns T \\0- *+(p, k) + i V \ \9 - $+{u, k) - irj\ ) 

For the factor \9 — $+(p, k) + ir/| _1 we proceed as follows 

• if 1 < I < l then we apply \9 - ®+(p, k) + ir]\~ l < Ce l ; 

• if I = Iq then we use 

\9-<S>+{p,k)+iT 1 \- 1 < I??!" 1 < Ce l °. 

Analogously, 



sup 

keQinSinSf 



(p)(u) 



(8.3) 

We finally obtain 

(8.4) sup 

ksQinSinSj 



$+{u,k) -ir]\- 1 < Ce l 



1<1<1 . 



1 



<3?+(p, k) + irj\ \9 - <J?+(it, k) - irj\ 



< Ce 



i+i 



\<1,1<Iq. 



2(1 



RICCARDO ADAMI AND LASZLO ERDOS 



Now we estimate su Pfe 6 Q. n 5 in g_ ( \p)(l) ) • Using inequalities (p!. I ). (3.2]) we have 



.5) (p)- 2 < C 



< C 



(kf 



(e(p + k)) ~ (6-<S> + (p,k)-9 + uj(k)) 



C 



(k)\6-<S> + (p,k)) 
(0) 



Notice that in Si, 1 < Z < Zo ° ne has (# — /c)) < 1 + e 'p, therefore such a quantity can be 

estimated by a constant and we end up with 



(8.6) 

Analogously, 
(8.7) 



(k) 2 



(u)- 1 < C 



keS u l < Z < Z . 



keSj,l<l<l . 



Therefore, from ([83D, §EM, (E2D one obtains ([82]) for 1,1 / 0. 

To estimate ZTj^o we observe that, proceeding like in (18. 5f) and using 



(6-$+(p,k)) 
' - $ + (p, k) +ir]\ 



< C, keS , 



one gets 
(8.8) 

Analogously, 
(8.9) 



®+(p,k) +irj\ (p) 



< C 



< c 



(k) 2 

(k) 2 



k £ S . 



k e S . 



Then, from (jEU) and one obtains i$EM for 1,1 = 0- Finally, to obtain i$872) for I = 0, I ^ it is 
sufficient to use inequalities (|8.7p . (|8.8p and (|8.3p . and analogously for the case Z / 0, Z = 0. 
The volume factor |Qj fl Si D 5r| in (18. ip is estimated as follows: 

• if both Z and Z are greater than zero, then we use inequality (13. 7p and estimate the volume of 
integration by 

„-l-l 



QiDSiD S T 



< C- 



\p — u\ ' 

if Z = and Z 7^ 0, then we use inequality (|3.6p and obtain 



n So n 



vice versa, if I > 0, Z = 0, then 



Qi n 5/ n 5 



If both Z and Z equal zero, then the volume 



< Ce~ l ; 



is estimated by Ce 



QiDSoD S 



is estimated by C. 
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Then, 



EE 



C(k) dk 



1 1,1=0 1 1 



< C 



- $> + (p,k) +irj\ |0- <$>+{u,k) -irj\(p)(u) 
I + 2l + l 2 /\p-u\^. ...... ^ n (log*r,) 2 



(0)2(0)5 Y (0)2(0)2 min(l,|p-u|) 

The decay property (|3.8p of £ guarantee that the sum in i is finite. Then, 

C{k)dk (log*r?) 2 



£ll<-> 4 £| 



< c- 



•10) 

" |0 - $+(p,fc) + irf\ \0 ~ $+(u,k) - ir]\(p)(u) (0)2(0)2 min(l, |p - u|) 

Inequality (|8.10p can be improved in the case \p — u\ < r/. Indeed, in this case, we can estimate terms 
with 1,1 > by means of inequality (|3.6[) in spite of (|3.7|) , We obtain 

£{k)dk 



.11 



|0 - ® + (p, k) + ir}\ |0 - $+(u, k) - irj\{p){u) 

log* 77 



<c 1 + 2/p + Zpe'Q ^ 

(0)2(0)1 Y 



> ^IIl°°(q,) < c 



(0)2(0)277 



Combining (f8T0|) and (jgTTl) we arrive at (j3~TT1) . 
Estimate p,12p is trivial. 



□ 
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